* Figure 2: INTERVENTION STOVE PURCHASE RATES

use "${output}panel_r0_r1_r2.dta", clear // Load merged analysis panel

********************************************************************************

* Keep households surveyed in each round
keep if sample_household == 1

* Label NGO and non-NGO households
label define ngo_status 0 "Non-NGO" 1 "NGO"
label values chirag_strata ngo_status

* Purchase rates (%) graph
clonevar purchased_intervention_stove_pc = purchased_intervention_stove
replace purchased_intervention_stove_pc = purchased_intervention_stove_pc * 100
cibar purchased_intervention_stove_pc if surveyround == 1 & treatment == 1, over(chirag_strata) ///
	level(95) ///
	bargap(10) ciopts(msize(0) lcolor(black)) ///
	baropts(base(0) yla(0(10)60) lcolor(black)) ///
	graphopts(ytitle("Percent of households") ///
	name(g_purchase, replace))

graph export "${results}figure_purchase_rates.png", width(5000) replace

gr drop _all
